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In this paper, we study a well known asymptotic limit in which the sec- 
ond Painleve equation (Pn) becomes the first Painleve equation (Pi). The 



■ limit preserves the Painleve property (i.e. that all movable singularities of 

o 

all solutions are poles). Indeed it has been commonly accepted that the 
\Q ■ movable simple poles of opposite residue of the generic solution of Pn must 

coalesce in the limit to become movable double poles of the solutions of Pi, 
even though the limit naively carried out on the Laurent expansion of any 
solution of Pi i makes no sense. Here we show rigorously that a coalescence 

o 

^ ' of poles occurs. Moreover we show that locally all analytic solutions of Pi 



arise as limits of solutions of Pn. 

1 Introduction 

An ordinary differential equation is said to be of Painleve type (or to possess the 
Painleve property) if the only movable singularities of its solutions are poles. The 
property is strongly related to integrable systems (systems which can be solved via 
related linear problems) [[l], |], IE EDI- Knowledge of equations with the Painleve 



property, including various methods of classification, is therefore valuable in the 
search for integrable systems. Asymptotic limits of differential equations that 
preserve the Painleve property provide another mechanism for such searches. 



In a series of papers published around the turn of the century, Painleve |TT[ . 
Gambier |J, and Fuchs || conducted an exhaustive search for all equations of 
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Painleve type of the form 

u" = $(x;u,u'), 

where $ is analytic in x and rational in u and v! . They discovered six equations 
of Painleve type which, up to a transformation, are the only equations of this 
type whose general solutions are new transcendental functions. These equations 
are known as the Painleve equations Pj-Pyi- The first two are 

Pi : u"{x) = 6u 2 + x; 

Pi i : u"(x) = 2u 3 + xu + a; 

where a is a complex constant. 

Painleve [12| noted that under the transformation 



x = e z — 6e~ . 



u = ey + e 
a = 4e~ 15 , 



2„ « £ -io. 

5. 



Pi i becomes 

y"(z) = 6y 2 + z + e 6 {2y 3 + zy] . (2) 

If e vanishes, equation (0), which we will refer to as Pn(e), becomes Pi with 
x replaced by z and u replaced by y. We will say that under flU), Pn degenerates 
to Pj and write Pn — > Pj. Painleve gave a series of such degeneracies which is 
summarized in Figure 1. 

Pin 
/ \ 
-Py/ — ► -Py -Pr/ — >■ Pi 

\ / 
P/y 

Figure 1: Degeneracies among the Painleve equations. 
Transformations which lead to these degeneracies can be motivated from the 



point of view of isomonodromy problems [?, [H|. Using a method based on max- 
imal dominant balances, Joshi and Kruskal J7| have found a new degeneracy of 
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Piv to another equation of Painleve type (equation XXXIV on p. 340 of Ince ||). 
Their paper raises the possibility of using asymptotic limits between differential 
equations which preserve the Painleve property as new tools in the search for, 
and classification of, integrable systems. 

The central concern of the present paper is an exploration of the convergence of 
solutions of -P//(e) to solutions of Pi as e vanishes. In particular we are concerned 
with the way in which simple poles of oppositely signed residues in solutions of 
Pn(e) coalesce to form the double poles of Pi. Unfortunately, a purely local 
analysis of this coalescence is problematic as the radius of convergence of any 
Laurent expansion centred on any pole necessarily decreases to zero if the poles 
coalesce. Rather than attempt to find an accurate upper bound on the radius of 
convergence we use techniques based on steepest ascent curves similar to those 
first expounded by Joshi and Kruskal in their direct proof that Pi to Pyi possess 
the Painleve property || . Before embarking on this problem we analyse a model 
problem given by a similar degeneracy between the autonomous versions of Pi 
and Pn whose general solutions are expressible in terms of elliptic functions. We 
show that the poles here coalesce by estimating the distance between them. Such 
estimates are obtained in section 2. 

In section 3 we consider general equations of the form 

-j 1 = Mz,Vi,---,Vn,e), l<i<n (3) 

where the fi are entire functions of (z, yi, . . . , y n \ e)- We will say that equations 
(|3D degenerate to the equations 

dy% 

-T = fi(z,yi,...,y n ;0), l<i<n (4) 

in the limit as e approaches zero. We show that, locally, any analytic solution of 
the target equations (^) can be obtained in the limit as e — > of a solution to 
equations (||). A corollary of the theorem states that if yi is a solution of Pj then 
given any compact subset K on which yi is analytic then there is a solution, y of 
Pn(e) such that y — > y t on K with respect to the sup norm as e — > 0. Hence, by 
considering the maximal analytic extension of y we see that y —>■ yi everywhere. 
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In section 4 we examine the rate of coalescence of poles. We obtain estimates of 
the distances between coalescing poles and show that these are of order e 3 . 



2 Two Autonomous Painleve Equations 

Consider the following autonomous versions of Pj and Pjj, 

Ej u" = Qu 2 + X 

u" = 2u 3 + fiU + a 

where A, fi £ C are constants and the primes denote differentiation with respect 
to x. The solutions of Ej and Ejj are either constants or may be expressed in 
terms of elliptic integrals. 

Following the analogy of the Pjj — > Pi coalescence we transform the variables 
in Ejj as follows: 

x = e 2 z, fi = Ae 2 - 6e~ 10 , u = ey + e" 5 , a = 4e" 15 . 

Under this transformation Ejj becomes 

jj = 6 y 2 + A + e 6 (2y 3 + Xy) , (5) 



where a dot denotes differentiation with respect to z, giving us the degeneracy 
Ejj — > Ej. In order to examine the nonconstant solutions of (|5]) we multiply the 
equation through by y and integrate. In this way we obtain 

f = e 6 P e (y) :=h + 2Xy + e 6 Xy 2 + Ay 3 + eV (6) 

where h £ C is a constant of integration. Take h given and fixed in the following 
analysis. The nonconstant solutions of equation (^) satisfy 

e 3 — = Q e (y) :-- 



d y ^ y/P^) 
Now, for e 7^ 0, 

P e (y) =■ (y - ao)(y - a>i)(y - a 2 )(y - a 3 ), 
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where 

«o = a (e) = -± + ^ + 0(e 12 ), (7) 
ai = ai(e) = r]i + 0(e 6 ), i = 1,2,3 

are the zeros of P e {y) and the r\i are zeros of Po{ff). 

We briefly recall some of the standard results from the theory of elliptic in- 
tegrals, beginning with a description of a Riemann surface for Q e (y) (see, for 
example Siegel ||14|| ). We will assume that h is such that for small e, P e {y) has 
distinct zeros (note that this is the generic case). Cut two nonintersecting slits 
in the Riemann sphere, say one from ao to a% and the other from a 2 to 03. Make 
two copies of the resulting manifold and label them OJti and WI2] these two slit 
spheres correspond to the two branches of the square root operation in the defini- 
tion of Q t {y)- Now take each side of both slits on 971 1 and identify them with the 
opposite sides of the corresponding slits of OK2. The resulting Riemann surface, 
9^, is homeomorphic to the 2-torus T 2 . Q e (y) is meromorphic throughout and 
the elliptic integral 

1(7) := / Qe(y)dy 

is well defined for any piecewise smooth curve 7 in D\ where y varies over the 
natural projection of 7 to the Riemann sphere CP 1 . 

Suppose that y has poles with residues of opposite sign at z + and z_, then 

z + -z_ = e- 3 [ Q € (y)dy, (8) 

where 7 is a path connecting 00 1 and oo 2 — the subscripts distinguish the points 
at infinity on the two slit spheres Wl\ and respectively. Such a path must 
pass through one of the open slits connecting the two spheres. Its projection onto 
the Riemann sphere must loop around the points a^, k — 0, . . . , 3 an odd number 
of times (note that if it encloses an even number of the points a^, the resultant 
integral is just a period of the elliptic function y). For small e the point ao is closest 
to infinity so we consider a path which begins at 00 1 and remains in DJli until it 
reaches the point ao, loops around it, and then retraces the corresponding path 
in WI2, terminating at 002- Since an arbitrarily small loop around ao contributes 
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nothing to (H), the distance between the two poles is simply 



\Z:—Z\ =2 



e- 3 



an 



Qe{y)dy 



(9) 



Next we refine our choice of the path of integration for the right hand side of 
equation (^|). At any point where y is analytic and neither y nor y' vanishes, there 
is a unique direction of fastest increase in \y\. Hence we can define a steepest 
ascent curve through any such point. A simple calculation using the Cauchy- 
Riemann equations shows that on such a curve, d\y\ = \dy\ (on a path of steepest 
descent, d\y\ = —\dy\). 

Let R + be the connected component of the region 

Q. := {z : \y(z)\ > \a \} 

containing z + in its closure. Note that the only pole in the closure of R + is z + 
(since there is a unique level curve of \y\ passing through every point in R+). 
Expanding y about a point z\ in the boundary of R + such that y(z\) = ao gives 

Viz) = a + -y"(z!)(z - z x f + 0((z- Zl ) 3 ) , 

where a ~ — 4e~ 6 and, from equation (||), y"(zx) ~ — 32e -6 . We see that z 1 is a 
(complex) saddle point. This implies that Z\ is the initial point for two steepest 
ascent curves (and two steepest descent curves). One of the steepest ascent curves 
must enter R + and terminate at z + . This is the path, T, over which we integrate 
in equation (|^). 

Choose r > so small that for all e such that |e| < r, 

|cio | > 2 max UaA} . 

l<i<3 

Note that this can be achieved because the expansion ([?[) shows that ao is large 
for small e. Then, since \y(z)\ > |oo| on T, for |e| < r and 1 < i < 3, we get 

\y\ < \y - Oil + \di\ < \y - Oil + \y\/2 \y\/2 < \y - a { \. 

So 

\ Q &)\ Z I == : „ • ( 10 ) 

y\v\ (\y\ - l a ol) 
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In particular, notice that the integral in (0) is convergent at infinity. 
Using (0), we find from equation ([]) that 

M 2V2 d\y\ 



\ z + 



< 



■2 k 



8V2 



l a o||N 3 \ 



y\ 3 (\y\ - l°o|) 

y=ao 



1 



|Qp| 

I y\ 



8V2 



y=oo 



l a o(ej 



0(e 3 



Therefore the two oppositely signed poles coalesce as e vanishes. 

In the above analysis we have only considered the generic case in which P e (y) 
has four distinct zeros. In the nongeneric case the Riemann surface, 9^, of y is no 
longer a torus. Our analysis, however, does not depend critically on the global 
topology of !H and the same estimates apply. 



3 Local Analytic Solutions 

The aim of this section is to prove the following theorem:- 

Theorem 1 Let (rji . . . , rj n ) be a given solution of the system of ODEs in 
which is analytic in some pathwise connected region Q C C and choose z e Q. 
Given any simply connected compact subspace K C Q containing zq, there exists 
a solution (y l7 . . . ,y n ) of equations (fjj and a number rx > such that, 

1. the yi are analytic in (z, e) for z £ K , |e| < r K ; 

2. yi (z,0)=r)i{z) WzeK; 

3. yi(z ,e) = rji(zo) Ve such that |e| < r^. 

Note that, regardless of the choice of K, the yi satisfy the same initial value 
problem at z . This theorem shows us that, locally, solutions of equations (|3]) 
converge onto solutions of equations @. It shows that the singularities of this 
family of solutions of equations (||) lie arbitrarily close to those of equations @ 
(or go to infinity), for small e. In the proof of this theorem given below we will 
make use of the following lemma which can be proved using elementary arguments 
involving majorant series (see, for example, Cartan 0). 
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Lemma 2 Consider the system of ODEs 

-p = fi(z,yi, ...,y„;e), l<i<n (11) 

together with the initial conditions 

y i (z 0> e) = <j)i(e), l<i<n 

where the (pi are analytic for |e| < r and i/ie /j are analytic on 

S := {0, . . . , y n ; e) : \z - z \ < p, \yi - (f>i\ < R, \e\ < r, l<i< n}. 

Then there is a unique solution y := (yi, . . . , y n ) of §T\ ) which is analytic in (z, e) 
whenever \e\ < r and 

\z-z \< Z^ R (e) := p ^1 — exp 

where 

M>sup|/i|, l<i<n. 
s 

Proof of Theorem 1: Since the fa are entire, we may expand them as power 
series, 

fi(z, yi,..., y n ; e) = E a )^-k n i^y k x ' ' ' 

which converge everywhere. 

Fix p,R,r Q > 0. Let T be any finite length curve connecting z to dK. We 
will first prove existence in a thin neighbourhood of T. Define 

B := sup{|z| : \z — z\ — p, z G T}, 
L := sup \i]i(z)\, 

l<i<n 

and M := maxKKn Mj where 

M t := 2Y,\a) kl ... k jB^R + L) k ^W . (12) 
This last series converges because the /; are entire. 
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'r)R 



(n + l)pM 



Let S = {(z,y 1 ,...,y n \e) : \z-z \ < p, |e| < r , ||/< — »fe(^b)| < 1 < 2 < n}. 
Then for (z,y 1: . . . ,y n ; e) G So we have 

|/i(2,yi,...2/ n ;e)| 

^ E \ a %u-kJ\ Z H\yi - mi*)] + \Vl(Zo)\) kl - ■ ■ (\Vn - Vn(z )\ + \Vn(z )\) 



kn I ^ 1 1 



< J2\a% u ,_, k jB^R + L) k ^r l 
= -Mi. 



Therefore sup zGSo \fi\ < M and so we deduce from Lemma |2| that there is a 
solution y(°) := (y^, . . . , y^) of equations ([3]) satisfying the initial condition 



y<°W)=iK*>), 



e < r . 



Furthermore, is analytic in (z, e) provided \z — Zq\ < Z^ ro R (e) (see Lemma | 
Notice that R (e) has the maximal value 

R \\ 



Let z\ be the first point on T such that \z± — zq\ = d/2 (if no such point exists 
then we have finished). 

Next we show that by restricting the range of e we can ensure that the initial 
value problem at z% gives us a solution whose radius of convergence in z is again 
bounded below by d/2. At z — zi, y^ -* is analytic in e for |e| < f for some f < ro- 
Let Si(f) := {(0,2/i,...,y n ;e) : 1 2; — 2r 1 < p, |e| < f, \yi- i q i {z Q )\ < R, l<i<n}. 
Then 

/ \ fei 



sup I fi I < 51 Kfcx-.fcJ^ LR+ sup 

Si(f) V 



e <r 



(0)/ 



R + sup 

v \e\<? 



(13) 



Now as r — > 0, supi e i <? \y\ (zi,e)| — ► |?7i(zi)| < L, and so ( fL3|) approaches ^M. 
Therefore there exists ri such that < r\ < f < ro and 



sup < M. 
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Invoking Lemma |2] again we see that there is a solution, y W , of equations (|3]) 
satisfying 

A,£) = y w (ft l £) 

for all |e| < ri, which is analytic in (z, e) provided \z — Zi| < Zff r R (e). We then 
look for the next point, z 2 , on T such that |#2 — z\\ = d/2 (if such a point exists) 
and repeat the above argument for a finite number of points z 2 , Z3, . . . , Zn in order 
to cover the curve. y( 4+1 ) analytically continues y^. y(z) := (yi(z), . . . ,y n (z)) is 
then defined to be whenever z lies in the domain of analyticity of y^. Since 
we proceed in steps of d/2 in z, the radius of convergence of y(z) about any point 
of r is bounded below. The compactness and pathwise connectedness of K then 
ensure that we can analytically extend y(z) to all of K by using a finite number 
of curves Tj from zq. The existence of the number then follows because we 
require only a finite number of reductions of r in the above analytic continuation 
of y(z). 

□ 

The condition that K be simply connected is essential for the single- valuedness 
of y(z). For example, consider the equation 

y = by + ez . 

The solutions of this equation for e = are elliptic functions and therefore mero- 
morphic. However, Painleve analysis (see ||) reveals that generic solutions to this 
equation for e 7^ possess logarithmic singularities. So locally analytic solutions 
of the equation with e = whose domain of analyticity (Q, in Theorem 1) is not 
simply connected do not necessarily arise from analytic solutions of the general 
equation, but rather from multivalued ones. 

In the case of Pi and Pn, however, the solutions are meromorphic Rill, Hi- So, 
on recalling the form of transformation ([!]), we see that all solutions of Pn{e) are 
meromorphic and therefore single valued. Analytically extending any solution of 
P//(e) along any path connecting zq to any other point will give a result which is 
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independent of the particular path chosen. Hence, when we apply Theorem 1 to 
Pn(e) we can weaken the requirement that K be simply connected, instead de- 
manding only that it be pathwise connected. The theorem then has the following 
corollary. 

Corollary 3 Choose zo,a,(3 G C. Let yj and y be maximally extended solutions 
of Pi and Pn(e) respectively, both satisfying the initial value problem given by 

y(z ) = a, y'(z ) = (3. 

Let Q C C be the domain of analyticity of yj. Given any compact K C ft, 
3rx > such that y is analytic in (z,e) for z G K, \e\ < and y — > yi with 
respect to the sup norm as e — > 0. 

Proof: Apply Theorem 1 using some compact pathwise connected subspace 
Kcil such that {z } U K C K. 

□ 

4 Coalescence of Poles and the Second Painleve 
Equation 

We now return to the problem of estimating the rate of coalescence of poles in a 
family of solutions to Pn{e) 

y" = 2eV + 6y 2 + e 6 zy + z, 

as e — > 0. Choose z G C. We will consider a family of solutions to Pn(e) given 
by y(zo) = a, y'(zo) = f3. Multiplying Pn(e) through by y' and integrating along 
some path 7 from z to z gives 

[y\z)f = eV + % 3 + 2zy + e 6 zy 2 - f [ly + e% 2 } dz + k t =: e 6 F e {z, y} (14) 

7 

where 

k e = I3 2 - {e 6 a 4 + 4a 3 + 2z a + e 6 z a 2 } . 
11 



From the corollary to Theorem 1 in section 3 we see that as e — > the 
solution to Pn{e) given by y(z ) = a, y'(z ) = (3 converges to the solution yj of 
Pi satisfying the same initial conditions, on any compact subset K of the domain 
of analyticity of yi. 

Suppose yi has a double pole at z. Let D be the closed disc of radius p centred 
at containing both z an d z in its interior. Let K be D after we have deleted 
open discs of small radius 5 centred at each of the poles of yi which lie in D. 
From Corollary 3 of Section 3 we see that for sufficiently small e, any simple pole 
of y which lies in D must be within 8 of a double pole of yi (since it cannot lie 
in K). 

Let z+,z~ G E>s(z) be the positions of two poles of y of oppositely signed 
residues. The distance between these poles is given by 

\z + — Z- 

for some path V between points whose natural projection to CP 1 is y — oo. The 
final integral in ([L5|) makes sense if we use the fact that locally on T, z may be given 
as a function of y. Also, since all solutions of Pn(e) are nonconstant meromorphic 
functions, the points on T at which F t {z, y} vanishes do not accumulate. 

As was the case in the coalescence of poles induced by the En — > Ei degen- 
eracy, the opposite signs of the residues of the poles of y at z_ and z + indicates 
that T must loop around a zero, z\ y of y' . The close proximity of the poles for e 
small indicates that y must be large at this stationary point. 

We take Y to loop around a zero, z\, of F 6 {z, y}, to be specified below. Let 
A := | y(zi) | . Define the region 

VI := {z G D : \y(z)\ > A}. 

Then Q is a union of regions surrounding poles of y. We take z\ so that the 
connected component, R + , of Q whose closure contains z + contains no other 
stationary points of y; i.e. y'(z) ^ for all z G R+. An analogous argument to 
that outlined in Section 2 shows that R + contains no pole other than z + and that 
Z\ is the initial point for two curves of steepest ascent and two curves of steepest 
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-7 r 


/ dz 




J 2_ 







dy 



Fe{z,y} 



(15) 



descent (each separated by one of the four level curves of \y\ which pass through 
zi). This follows from the fact that y(zi) ^ 0, y'{z\) = 0, and y"(zi) ^ 0. One 
of these steepest ascent curves, r + say, lies in R + and so connects z\ to z + . The 
other steepest ascent curve, T_, lies in Q \ R + and is of finite length (necessarily 
terminating at a pole, z_ say) since, for large A, Q is a union of small disjoint 
regions containing the poles of yi. We take the path of integration, T, in equation 
( |15D to be the union of these two paths. 

Since z\ is in the boundary of Q and is the initial point for a curve of steepest 
descent, there is a curve connecting z to Z\ contained in D \ £1. The (initial) 



path of integration, 7, in equation (|T4]) connecting zq to z G T is taken to be this 
descent curve followed by one of the steepest ascent curves, Ti or T 2 , from z\ to 
the point z. 

We now estimate F € {z, y} for large y on such a curve. For z G T, we have 

sup |s/(OI = \y( z )\- 

Since D\Q only contains a finite number of small holes, the length of the path 7 
from zo to any point on T can be bounded by some d > which is independent of 
e for e small. Also, for small e, k £ can be bounded above by c 2 , say, where c > 
is independent of e. 

From equation flT^|) we see that A := |y(2i)| is asymptotically close to 4|e|~ 6 . 
Let r > be an upper bound on e 6 which is so small that 

r < max |c? _1 , P -1 } , and A > max {c, <i, p} . 

We now see that on T, for |e| 6 < r, 

e 6 F e {^} = e y + % 3 + <K^), 



where 



< \2zy\ + 



e zy 











+ 


J Z 2ydz 


+ 






7 




7 



+ I fee 



< 2p|y| +rp|?/| 2 + 2d|y| +rc%| 2 + c 2 
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where k = 5 + r(d + p). 

So on T, 4>(z,y) = ip(z,y)y 2 where 



< EflVlIlE 



e b F e {z,y} = eY + W + ^(z,y)y- 

Now arg(y) is a constant along any path of steepest ascent for y 
d\y\ = \dy\ there). Hence r + can be parameterized by t £ (1, oo), where y 
yi := y(z\). Since F e {z\,yi) = 0, we see from ( |16D that 



(16) 

^since 
= tyi, 



+ 1 



as e — > 0, ffivinc 



lim 



e G yi 



,0 4 



-1. 



(17) 



So if we hold t fixed as e — > we have 



lim — - = —t. 



Consider the following ratio 

F e {z,y} 





2 

) + 


m 


e 6 ib 

+ 4 




2 


, 4 , 







Using the limits in fll?] ) and (18), we see that 



limi?, = 1. 

The definition of limit then shows for some given < u < 1, 3r>0 sufficiently 
small such that for all e with e 6 < r, 

\F e {z,y}\ > v 2 \y 3 (y- yi 

Since the same argument holds on T_, we have from equation (|T5|) 



\Z-l. — Z- 1 < 



2/1 



/o ° \/y 3 (y-yi) 

where the integration is along a path of steepest descent (from a pole of y). So, 
recalling that along such a path, \dy\ — —d\y\, we have 

i 2 f A d \y\ 4 

\z + — Z-\ < ' 



?/ !'T' J \y\ 3 (\y\ - A) vA \ e \ 
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Since A ks 4|e| 6 , this shows us that the distance between the poles of solutions 
of Pn(e) is of order e 3 . 
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